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Abstract 



The effect of multi-particle Coulomb final state interactions on higher-order intensity 
correlations is determined in general, based on a scattering wave function which 
(^ ■ is solution of the n-body Coulomb Schrodinger equation in (a large part of) the 

^ , asymptotic region of n-body configuration space. In particular, we study Coulomb 

effects on the n-particle Bose-Einstein correlation functions of similarly charged 
0> ' particles and remove a systematic error as big as 100 % from higher-order multi- 

particle Bose-Einstein correlation functions. 
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1 Introduction 



One of the most fundamental quests of high-energy physics is the determi- 
nation of the phase diagram of strongly interacting matter. At high densities 
and/or temperatures the quarks are expected to be liberated from their con- 
finement within hadrons and a new phase of matter, the Quark Gluon Plasma 
(QGP), to be formed. In order to explore this new phase of matter, the Rela- 
tivistic Heavy Ion Collider (RHIC) has been constructed at Brookhaven Na- 
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tional Laboratory, to collide Au + Au nuclei at ^/s = 200 AGeV center of 
mass energy. 

One of the new features of RHIC physics will be the production of 600 - 1200 
charged pions per unit rapidity. Due to this reason, the PHENIX and STAR 
detectors will be able to determine not only single particle spectra and the two- 
particle Bose-Einstein correlations, but also the higher-order Bose-Einstein 
correlation functions, which turn out to be essential in order to distinguish 
fully chaotic and partially coherent particle sources from each other [1]. As 
partial coherence is a fundamental aspect of quantum fields, and it can be 
related to a possible Bose-Einstein condensation of pion wave-packets to the 
wave-packet with the smallest energy in the rest frame of the source [2], the 
determination of the higher-order Bose-Einstein correlation functions at RHIC 
is of great theoretical interest [1,3-9], as well as a great experimental challenge. 

However, Coulomb (and possibly strong) final state interactions of pions play 
an important role in shaping the final multi-particle Bose-Einstein correla- 
tion functions. As no consistent and systematic treatment of the final state 
interaction of a charged multi-boson system is available in the literature, the 
experimental removal of the Coulomb effects from the n-particle Bose-Einstein 
correlation functions is based hitherto only on some ad-hoc generalization of 
the Gamow formula to the multi-particle case. 

In the present Letter we propose a straight-forward method for a systematic 
quantum-mechanical treatment of Coulomb final state interactions in higher- 
order Bose-Einstein correlation functions. Although the validity of the method, 
described below, is limited to a certain, albeit large kinematic domain {Q,q ), 
due to the fact that the exact solution of even the 3-body Coulomb scattering 
problem is beyond presently available means, we think that the results pre- 
sented here represent a first important step towards establishing a link between 
few-body physics and Bose-Einstein correlations in high-energy multi-particle 
physics. Especially, the new Coulomb wave-function corrections indicate that 
the generalized Gamow correction method would make a factor of two error 
in the 5-th order Bose-Einstein correlation functions, if the radius parame- 
ters were in the 5 - 10 fm range, as characteristic for heavy ion collisions. If 
the characteristic radius parameters were as small as 1 fm, the characteris- 
tic size in reactions of high-energy particle physics, the generalized Gamow 
factors would be acceptable, at a 10 % level of precision in the 5-th order 
Bose-Einstein correlation functions. 

The strength of higher-order correlation functions increases much slower for 
partially coherent particle sources than for incoherent sources with an unre- 
solvable halo of long-lived resonances [3,1]. Even if the strength of the second- 
and third-order Bose-Einstein correlations were similar in a partially coherent 
and another incoherent particle source, the strength of the 5-th order correla- 



tion functions would be a factor of 2 different between the partially coherent 
and the fully chaotic cases [1]. In order to distinguish these scenarios, the 
Coulomb final state interactions must be corrected for, and the error on the 
Coulomb correction must be kept under control. 

In a recent Letter, we have presented a refined treatment of the 3-body 
Coulomb correction problem [11], with application to new high-energy heavy- 
ion data by the NA44 experiment [12]. In this work, we generalize this 3-body 
Coulomb wave function integration method to the case of n-particle Coulomb 
corrections. 



2 Bose-Einstein n-particle correlations and final state interactions 

Let us summarize some basic properties of the Bose-Einstein n-particle corre- 
lation functions using only the generic aspects of their derivation, and establish 
a link between the theory of final state interactions in few body physics and 
the theory of Bose-Einstein correlations in high-energy particle and nuclear 
physics. 

The n-particle Bose-Einstein correlation function is defined as 
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where iV„(ki, ■ ■ ■ , k„) is the n-particle inclusive invariant momentum distribu- 
tion, while A^i(ki) is the single-particle invariant momentum distribution. It is 
quite remarkable that this complicated object, that carries quantum mechan- 
ical information on the phase-space distribution of particle production as well 
as on possible partial coherence of the source, can be expressed in a relatively 
simple, straight-forward manner both in the analytically solvable pion-laser 
model of Refs. [13,2] as well as in the generic boosted-current formalism of 
Gyulassy and Padula [14] as 
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where o"*^"^ stands for the set of permutations of indices (1, 2, ■ ■ ■ , n) and ai 
denotes the element replacing element i in a given permutation from the set 



of (T^") , and, regardless of the details of the two different derivations, 

G'(k„k,) = (at(k,)a(k,)) (3) 



stands for the expectation value of a^(ki)a(kj). In the boosted-current for- 
malism, the derivation is based on the assumptions that (i) the bosons are 
emitted from a semi-classical source, where currents are strong enough so 
that the recoils due to radiation can be neglected, (ii) the particle sources 
are an incoherent random ensemble of such currents, described by a boost- 
invariant formulation [14], and (iii) that the particles propagate as free plane 
waves after production. However, a formally similar result is obtained when 
particle production happens in a correlated manner, and even final state inter- 
actions between the produced particles are allowed for, generalizing the results 
of Refs. [2,15,16]. 

In the pion-laser model, the n-particle exclusive invariant momentum distri- 
butions read as 

n 
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with 

G'i(k,k,)=Tr{/)iat(k,)a(k,)}, (5) 



where pi is the single-particle density matrix in the limit when higher-order 
Bose-Einstein correlations are negligible. One can show [15,9], that the n- 
particle inclusive spectrum has a similar structure, if the multiplicity distri- 
bution is Poissonian in the rare gas limit: 
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The functions G„(kj, k^) can be considered as representatives of order n sym- 
metrization effects in exclusive events, see refs. [13,2,15] for more detailed 
definitions. The function G(kj, kj) can be considered as the expectation value 
of a^(kj)a(kj) in an inclusive sample of events, and this building block includes 
all the higher order symmetrization effects. In the relativistic Wigner-function 
formalism, in the plane wave approximation G(ki,k2) can be rewritten as 



G(ki, k2) = / d^x S{x, K12) exp(igi2 ■ x) (8) 

i^i2 = 0.5(A;i + A;2) (9) 

qi2 = h-k2, (10) 



where a four- vector notation is introduced, k = (ym^ + k^, k), and a-b stands 
for the inner product of four-vectors. Due to the mass-shell constraints, i.e. 
Ey_ = \Jvn} + k^, G depends only on 6 independent momentum components. In 
any given frame, the boost-invariant decomposition of Eq. (8) can be rewritten 
into the following, seemingly not invariant form: 

G'(ki,k2)= / (i^x ^Kialx) exp(?qi2x), (11) 

'S'ki2(x)= / dt expiifB^.^qut) S{x,t,Ki2), (12) 

/3^^^ = (ki+k2)/(i?l + i?2). (13) 

Note that the relative source function Ski2 (x) reduces to a simple time integral 
over the source function S{x, K) in the frame where the mean momentum of 
the pair (hence the pair velocity /3^^J vanishes. 

If n particles are emitted with similar momenta, so that their n-particle Bose- 
Einstein correlation functions may be non-trivial, Eqs. (5,7) will form the 
basis for evaluation of the Coulomb and strong final state interaction effects 
on the observables. On this level, all the correlations are build up from cor- 
relations of pairs of particles. This is due to the specific form of the density 
matrix that includes just the right amount of stimulated emission to make 
a further calculation straight-forward. Note also that a similar result can be 
obtained in the semi-classical boosted-current formalism, where the particle 
production has negligible effect on the elementary source of pion production 
(bremsstrahlung-like radiation). 

Let us point out that the exact solution of multi-particle Bose-Einstein sym- 
metrization in the pion laser model resulted in a Poisson cluster picture [15,2]. 
This implies that in the rare gas limit, without Coulomb or other final state 
interactions, the multi-boson correlations appear only as random admixture 
of a small amount of correlated pairs to independently distributed single par- 
ticles. As the density increases, also the fraction of correlated pairs increases 
and the admixture of independently distributed clusters of particle triplets, 
quartets, and higher-order n-tuples becomes correspondingly more important. 
The result of refs. [15,2] indicates that below the onset of Bose-Einstein con- 
densation, a fully symmetrized multi-boson system can be considered as a 
convolution of independently distributed clusters of particle n-tuples, and it 
is natural to apply Coulomb corrections within such clusters of particles only. 
We shall also discuss that, when one of the particles becomes separated from 
its cluster, the relevant n-particle Coulomb correction factor will reduce to 



the Coulomb correction factor of a smaller cluster that contains the remaining 
n — 1 particles. 



3 Quantum mechanical treatment of the Coulomb n-body problem 



In order to treat correctly the Coulomb corrections to the n-particle correla- 
tion function, knowledge of the n-body Coulomb scattering wave function is 
required. We restrict ourselves to the case that the transverse momenta of all 
the particles in the final state in their center of mass are small enough to make 
a nonrelativistic approach sensible. Hence the problem consists in finding the 
solution of the n-charged particle Schrodinger equation when all n particles 
are in the continuum. 

Consider n distinguishable particles with masses rrii and charges Cj, i = 1,2, 
■ ■ ■, n. Let Xj and kj denote the coordinate and momentum (three-)vectors, 
respectively, of particle i. From these we construct in the usual manner the 
relative coordinate r^ = Xj — Xj and the relative momentum kj^ = (m^kj — 
mikj)/{mi + rrij) between particles i and j, the corresponding reduced mass 
being /i^ = mimj/{mi + rrij). 

The n-particle Schrodinger equation reads 



Ho+Y. y^^ - E \ <.!.k„(xi, ■ ■ ■ , x„) = 0, (14) 



where 
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is the total kinetic energy for n particles in the continuum. Hq is the free 
Hamilton operator and 



V^Jir^J) = Vi^(r,,) + Vg(v,,) (16) 

the interaction potential between particles i and j, consisting of a strong but 
short-range {V^) plus the long-range Coulomb interaction (V^j^(rjj) = e^ej/ \ 
Tij I). Equation (14) has to be complemented by the complete set of boundary 
conditions in order that a unique solution be obtained. 

Already for n = 3, the exact numerical solution of the Schrodinger equa- 
tion (14) for ii^ > is beyond present means, partly for principal and partly 



for practical reasons. For a brief discussion of the related difficulties see [17]. 
But, at least the complete set of boundary conditions to be imposed is nowa- 
days known analytically [18,19], in the form of the explicit solutions of the 
Schrodinger equation in all asymptotic regions of the three-particle configura- 
tion space. Apart from the trivial two-cluster region relevant for an asymptotic 
configuration containing only two particles one of which is a bound state of 
two particles, the asymptotic solution takes its simplest form in the asymp- 
totic region conventionally denoted by Qq and characterized by the fact that 
- roughly speaking - all three interparticle distances become uniformly large, 
i.e. all I Tij 1^ oo (for a precise definition of the various asymptotic regions 
see [19]). There exist three more asymptotic regions Qij, i < j = 1,2,3, which 
are pertinent to situations characterized by final state interactions between 
particles i and j. But the appropriate asymptotic solutions are rather more 
complicated. From the physical point of view the union of all these regions 
^0 U ^12 U ^13 U ^23 is relevant for the complete breakup into three free parti- 
cles. As has been shown in [11], in spite of the lack of an exact solution of the 
three-body Schrodinger equation in the whole three-body configuration space, 
already knowledge of the asymptotic solution in Qq led to a systematic, well 
controlled extraction of Coulomb effects in the three-particle Bose-Einstein 
correlation measurements, in contrast to earlier, ad-hoc 3-body Coulomb cor- 
rection methods. 

In the final states of heavy-ion reactions, where a large number of charged 
particle tracks appear, the mutual, macroscopically large separation of tracks is 
one of the criteria of a clean measurement. This suggests that in order to study 
n-body correlation functions, again knowledge of the wave function in fig , 
the region in n-particle configuration space where all interparticle distances 
become uniformly large, i.e., | rj^ \^ oo for all values of {ij), may be sufficient. 
Here, uniform divergence of interparticle distances in Hq means roughly that 
<\ Tij \ / \ Tki \< oo for asymptotically large times for any two arbitrarily 
chosen particle pairs, although the interparticle distances themselves diverge 
for any pair. One immediate consequence of this definition is that in Oq 
the short-range interaction parts V^'j play no role any longer and can thus be 
neglected. 

For want of an exact n-particle Coulomb scattering wave function an approx- 
imate solution of Eq. (14) is sought. For this purpose, let us introduce the 
continuum solution of the two-body Coulomb Schrodinger equation as 






i 2,., ■^«'-^"'- 2^1 <*'<-■■''-«• <"' 

describing the relative motion of the two particles i and j with energy k?/2/ijj. 
The explicit solution is 



'^i:,iJ^j)=N^J e'^>^"'^F[-^77.j,l;«(| k,,- || n, \ -^-r,,-)], 



with Nij = e~'^'^iil'^ r(l + ir]ij\ and riij = CiCj^ij/ \ rij \ being the appropriate 
Coulomb parameter. F[a,b;x] is the confluent hyp ergeo metric function and 
r(a;) the Gamma function. With it the following ansatz for an approximate 
n-particle Coulomb wave function is made 



<:...,k„(xi,..-,x.) ~ n ^kr(r..)- (19) 



This ansatz can be justified by the following arguments. 

(i) The wave function (19) is asymptotically correct in the asymptotic region 
Oq ; that is, it is the leading term if all interparticle separations go to in- 
finity, of the (unknown) exact solution of the Schrodinger equation (14) [18]. 
Of course, for non-asymptotic particle separations it represents a theoretically 
not compulsory though plausible extrapolation. 

(ii) In the formal, time- dependent scattering theory the basic object is the 
M0LLER operator which maps the free n-particle state onto the correspond- 
ing scattering state. The mathematically rigorous definition of the n-charged 
particle M0LLER operator [20] requires, in contrast to the case of purely 
short-range interactions between the particles, the introduction of a 'renormal- 
ization factor'. The latter has the form of a product of n{n — l)/2 renormal- 
ization factors each of which is appropriate for the definition of the M0LLER 
operator for one of the possible pairings of the charged particles. Obviously, 
the ansatz (19) of the n-particle wave function is consistent with this renor- 
malization prescription. 

(iii) By suitably decomposing the (stationary) M0LLER operator of an n- 
body system into a chain of M0LLER operators of subsystems with fewer 
interacting particles, a wave function of the type (19) has been suggested as a 
lowest-order term of an n-particle Coulomb wave function in [21] to be used in 
all of configuration space. The assumptions entering were neglect of genuine 
higher-than-two particle correlations in the wave function, which is justified in 
Qq , and restriction of all two-particle scatterings onto their respective energy 
shells. 

(iv) The wave function (19) coincides for any selected particle triplet with the 
form pertinent to the given pre-selected triplet, if the corresponding interpar- 
ticle distances diverge [11]. 

(v) For n = 3, such an approximate wave function has been proposed in 
Refs. [23,24]. Although it ceases to be solution of the Schrodinger equation 
(14) for non- asymptotic values of the relative coordinates, it is nevertheless 
widely used, with considerable success, to calculate cross sections for the ion- 
isation of hydrogen atoms by (energetic) electron impact. 



The foregoing discussion makes it clear that the ansatz (19) is justifiable only 
for sufficiently large interparticle separations, a condition which is not easily 
translated into an experimentally accessible criterion. However, for n = 3 such 
a criterion has been established, namely that the total kinetic energy EU^i of 
three particles of unit charge be at least 0.2hc/Rc{ii^) which equals 10 MeV 
for a typical source size Rg = 4 fm [11]. Thus, assuming also for an arbitrary 
number n of particles (again for simplicity taken to have unit charge) the 
total kinetic energy E^^^^i to be equally distributed over the relative kinetic 
energies between each pair, the latter criterion generalises to the condition 
-Ejotii — 0.033n(n — l)hc/ Rdiva.) MeV. Hence, for the following we assume 
the n-particle Coulomb wave function to be given everywhere as 

^S..k„(xi,---,x„) ^ ^Wkjxi,---,x„) 

n 

:= ^jW) n <.''V..), (20) 

for Eit > 0.066 ^^%~ ^^ „ ^' MeV, (21) 

where M^"^' is an undetermined overall normalization constant. This is the 
building block for a properly symmetrized n-body wave function where the 
bosonic or fermionic nature of any subset of identical particles has to be taken 
into account in the symmetrization (or anti-symmetrization) process in the 
standard manner [22]. 

In this Letter we present explicitly the fully symmetrized wave function only 
for the case of n identical charged bosons, as to our knowledge measurements 
in high-energy physics attempting to reveal the strength of the multi-particle 
Bose-Einstein correlation effects[12,10] exist only for this special case. 

The fully symmetrized n-particle wave function has the form 

^Lt^.'^Jxi, . . . , x„) = 4^ 5: ^Lt.^..k„(x.„ ■ ■ ■ , X.J, (22) 

where o"*^"-* stands for the set of permutations of n different indices, and o"j for 
the permuted value of the index i in one of the permutations belonging to the 
set o"*^"-*. Using the ansatz (20), the above equation simplifies as 



*t:.'i (X, . . . . . X,.) = ^ i: n <*' (-..,,), (23) 

which contains only the two-body relative Coulomb wave functions. 
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The physics of the above ansatz is very simple: if all n final charges emerge 
into the continuum and if all are well separated from the other tracks, only 
the pairwise Coulomb relative wave functions play a role. However, relative 
Coulomb wave functions have to be taken into account for all possible particle 
pairs as the Coulomb interaction is of long range. Graphically, if we represent 
the n particles by n crosses, the relative Coulomb wave function between 
particle i and j can be represented by a line connecting cross i with cross 
j, and the full, asymptotically correct n-particle Coulomb wave function is 
represented by connecting each of the n crosses with the n — 1 others by 
forming a polygon with n corners and n{n — l)/2 lines (diagonals and edges). 

One can apply a simple approximation to (23) which preserves at least some 
features of the Coulomb distortion effects. It consists in neglecting, for any two- 
particle Coulomb wave function ■j/'j^ (^jj), the hypergeometric function in the 
exact solution (18) and retaining only the part e^^'^'^^^ Nij . After evaluating the 
double sums over all permutations of o"„ in a product, one finds 
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*L%,.(x,, 



(24) 



In the last line, the symmetrized n-particle wave function for neutral par- 
ticles \1/|^^ ..i^^(xi, ■ ■ ■ , x„) with its own normalization constant A/q has been 
introduced. As usual. 



G„- := |A,,,f = e-^"- \r{l +tr],j)f 



'ij 



"iji 



(25) 



is the Gamov penetration factor for the particle pair (ij). 

It, thus, follows that the proper generalization of the Gamow penetration 
factor for n charged particles reads as 



Gi .... 






(26) 



This expression contains n{n — l)/2 factors, corresponding to all possible pair- 
ings (ij). Moreover, it is self-consistent: if the momenta k',k",- ■ ■ ,k^^^ of / 
particles approach infinity such that for no two momenta /;;'•*•' and k^^^ the 
corresponding relative momentum remains finite, we have 



lim 



lim Gi^ 



^ai,---,i 



(27) 
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where the remaining n — I particles whose momenta remain finite are denoted 
symbohcally by cti, ■ ■ ■ , «„_;. Specifically, 

lim Gi,.,, = G'i,..,(„_i). (28) 

kn—>00 



An explicit check for n = 4 shows that, indeed, 

lim 6*1,2,3,4 = lim G12G13G14G23G24GU = G12G13G23 = 61,2,3- (29) 

A;4— >oo A;4— >oo 



Hence, the generalization (26) of the Gamow correction factor to arbitrary 
values of n is done in a self-consistent manner that satisfies its physically 
expected reduction property. 

This result was substantiated for the case of n identical charged bosons. In 
general, the final state of a high-energy heavy-ion reaction contains many dif- 
ferent kind of particles, with different charges and quantum statistical proper- 
ties. Nevertheless, the ansatz given in Eq. (20) can be utilized for any values 
of the charges, and the result can be symmetrized for a generic mixture of 
particles as prescribed in Ref. [22] . 

Let us add two comments: 

(i) As mentioned above, a wave function of the type (19) implies that the 
relative motion of each of the pairs of particles is independent of that of the 
other pairs, i.e., that no correlations between the motions of the particle pairs 
occur. In other words, the proposed form of the factorized n-particle Coulomb 
wave function does not include genuine higher-order correlations, only those 
that can be built up from two-particle Coulomb correlations. This is the same 
level of approximation that is used to derive Eq. (2), the generic form of the n- 
particle Bose-Einstein correlation functions. Results of Ref. [2,15] suggest that 
Eq. (2) is valid only if the density of bosons is below the limit of Bose-Einstein 
condensation. 

(ii) It should be kept in mind that the extrapolation out of the region fig 
implied by (19) is highly non-unique. Even for the 3-body Coulomb problem, 
various different wave functions which, of course, coincide asymptotically in 
Qq with (19) have been, and are still being, developed. 



4 Application to high-energy heavy-ion and particle coUisions 



The correlation function measuring the enhanced probability for emission of n 
identical Bose particles is given by Eq. (1). This correlation function is usually. 



11 



due to meager statistics, only measured as a function of the Lorentz invariant 
Q„, defined by the relation 



Q^ = E ?: 



i<j=i 



(30) 



where qij = h — kj, and where ki is the four-momentum of particle i. 
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Fig. 1. Filled circles stand for the Coulomb correction factor of 5-particle Bose- 
Einstein correlation functions for a source size of i? = 1 fm as obtained from the 
numerical integration of the 5-body Coulomb wave function, while the squares indi- 
cate the results of the less substantiated 5-body Gamow corrections; the inset shows 
the ratio of these two correction factors. Lines are shown to guide the eye. 

We can now calculate the Coulomb effects on the n-particle correlation func- 
tion using 



K, 



Coulomby^n 



{Q. 



/niLirf'xip(xi) 



^K„(x. 



(0)5 



/n?=ic/3xip(xi) </.^,jxi,..-,x 
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(31) 



where p(xi) is the density distribution of the source for particle i, taken as a 
Gaussian distribution of width R in all three spatial directions. This formu- 
lation makes it possible to extract information on the source size i?, and to 
compare this value with that extracted by means of a generalized n-particle 
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Fig. 2. Same as Fig. 1 but for n = 4 particles and i? = 5 fm 



Gamow approximation through K}jJ^i^^^{Qn) = nr<j=i ^ij- To this purpose 
we use the NA44 data sample of three pion events produced in S-Pb colhsions 
at CERN [12]. 

We have calculated the Coulomb correction factor, i.e. -ft'co«iom6('5n) [12], for 
source radius values i? = 1, 5, and 10 fm, for n = 2, 3, 4, and 5 particle 
correlations. The radii were chosen to be in the range of interest for high- 
energy particle and high-energy heavy-ion physics. The results are compared 
to the generalized Gamow approximation. We have checked that in the limit 
i? ^ the n-particle Gamow approximation is indeed recovered numerically. 

In case of a characteristic 1 fm effective source size, typical for Bose-Einstein 
correlation functions in various elementary particle reactions, the difference 
between the n-particle Gamow and Coulomb wave function corrections were 
smaller than 10 % for n = 4 and 5 particles, the n = 5 case being shown in Fig. 
1. However, for future measurements of 5-particle Bose-Einstein correlations 
in particle physics that aim at a precision better than 5 percent relative error. 
Coulomb wave function integration will be a necessity. 

For source sizes of 5 or 10 fm, that are the characteristic expectations for Au 
-|- Au reactions at RHIC, the difference between the results of the Gamow and 
Coulomb wave function corrections increased dramatically, see Figures 2-4. 
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Fig. 3. Same as Fig. 1 but for n = 5 and i? = 5 fm 



We find that for a source radius of 5 fm, we need to take this detailed calcu- 
lation into account already for the precise determination of the three-particle 
correlation function. However, with increasing number of particles, the devia- 
tion between the n-particle Gamow and the n-particle Coulomb wave function 
integration method increases drastically. For 5-particle Coulomb correction, 
the better substantiated Coulomb wave function integration method yields a 
factor of 2 deviation from the naive generalized Gamow method. 

At the end we note that systematic improvements of our treatment are pos- 
sible by i) including also effects of strong interactions between the particles 
of each pair, ii) replacing the simple product of Gaussians by a more realis- 
tic model for the production of particle n-tuples, and iii) invoking improved 
n-body Coulomb wave functions that are correct in a larger region of n-body 
configuration space than 0,^ ■ 

However, corrections i) and ii) are estimated to be small, as 1) the final state 
Coulomb interaction dominates over the final state strong interaction due to 
its long range and the relatively large source size, and 2) the effective source 
of particles is known to be approximately Gaussian from detailed studies of 
two-pion correlations. Correction iii) is also expected to be small, as a clean 
measurement of particle n-tuples will likely require that these particles be in 
Hq ; however, one has to wait till 4-th and 5-th order correlations are measured 
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in heavy-ion collisions in order to determine the more detailed experimental 
conditions. 



5 Summary and conclusions 



On the basis of an explicit, analytically given form of the n-body Coulomb 
wave function that is - at least asymptotically - correct in a large region of n- 
body configuration space, we have developed a new method to systematically 
correct for explicit many-body Coulomb effects which is applicable to data 
analysis in a broad range of measurements in high-energy particle and heavy- 
ion physics. A generalized Gamow correction factor has been established as a 
limiting case of vanishing source sizes. 

Specifically, we have worked out our approach for 3, 4, and 5 identical charged 
particles and have tested it for Gaussian source sizes with i? = 1, 5, and 10 
fm. We have numerically found that the generalized Gamow approximation is 
not reliable enough to determine the magnitude on the 5 % level of the five- 
body Coulomb correction factor if i? = 1 fm, the characteristic lenght-scale 
of strong interactions in high-energy particle physics. The range of interest 
in high-energy heavy-ion physics was probed in the R = 5 and 10 fm cases. 



15 



and systematic errors, as large as 100 %, were shown to be generated with 
the earher Coulomb correction techniques for the correlation function of 5 
particles. 
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